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Abstract 

We discuss scenarios in which the galactic dark matter in spiral galaxies is 
described by a long range coherent field which settles in a stationary con- 
figuration that might account for the features of the galactic rotation curves. 
The simplest possibility is to consider scalar fields, so we discuss in particular, 
two mechanisms that would account for the settlement of the scalar field in a 
non-trivial configuration in the absence of a direct coupling of the field with 
ordinary matter: topological defects, and spontaneous scalarization. 
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I. INTRODUCTION 



It has been known for a long time that the motion of the stars and gases 
around the center of most galaxies can not be explained in terms of the lu- 
minous matter content of the galaxies, at least not within the context of 
Newtonian gravity (see Q for a review) . The standard view is that there is in 
almost every galaxy a large component of non-luminous matter (the gravita- 
tional dark matter) that forms an halo around the galaxy and that provides 
the additional gravitational attraction needed to explain the "rotation curves" 
in terms of standard gravitational theory. There are several proposals for this 
dark component, ranging from new exotic particles such as those predicted by 
supersymmetry ||2|-Q, to other less exotic candidates such as massive neutri- 
nos, all collectively known as WIMPs (weakly interacting massive particles) 
(see 0j5| for a review), to the relative mundane idea of dark but ordinary 
bodies such as Jupiter-like objects collectively known as MACHOs (Massive 
Compact Halo Objects) p. Searches for these types of objects have been 
made 0, and although they report some findings, there doesn't seem to be 
enough of these objects to account for the galactic dynamics. Moreover there 
are severe bounds on the amount of baryonic matter in the universe arising 
from big bang nucleosynthesis and for some values of the Hubble constant 
those bounds also imply that some of the galactic dark matter ought to be 
exotic ||8-10|. Independently of this and despite their popularity, these type of 



models suffer from various problems and require surprising coincidences (see 
for example [Q). 

Another type of proposal, which is in some sense more radical, is based on 
the idea that the gravitational theory would have to be modified when deal- 
ing with the scales associated with the motion of stars in galaxies |12,13|, in 



particular the idea is embodied by the proposal of Milgrom that the laws 
of motion are modified when the accelerations involved are extremely small. 
Unfortunately this scenario has not, as yet, been converted into a fully rela- 
tivistically invariant theory. Another type of model that has been exposed is 
to replace general relativity by a higher order in curvature theory, which in 



some particular cases appears to be obtaining encouraging results |15]. The 
problem with this approach is that these types of theories have in general 
problems of principle like for example the lack of a well posed initial value 
formulation. Nevertheless such relatively radical proposals are still attractive, 
due in part to an intrinsic problems of the more conservative approaches in 
explaining the generality and universality of the phenomena, namely the fact 
that the amount of luminous matter seems to be such a good indicator of the 
amount of the dark matter component and the fact that the dark com- 
ponent happens always to distribute itself in such a way that the resulting 
rotation curves (hereafter referred to as RC) are almost flat away from the 
galactic centers [p 
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Thus, in contrast with the former scenarios which would need to assume 
not only the existence of the dark matter but also give some evolutionary 
scenarios that result in the aforementioned universality in its distribution, 
the modified gravity scenario would naturally account for such correlations 
without the need for additional assumptions. On the other hand the former 
scenarios do not present any problem in lending themselves to an acceptable 
theoretical formulation, compatible with present theories of particle physics 
and general relativity. 



The object of this article is to discuss a third type of scenario which has 
some of the advantages of each scheme. The idea is to take dark matter to 
be described not by a bunch of particles whose distribution needs to be ex- 
plained but by a coherent field which would settle in a universal stationary 
configuration that would account for the generic features of the RC. The sim- 
plest possibility is provided by scalar fields, which would of course have to 
be very long ranged (i.e. masses smaller than where Rq is the radius of 
the largest galaxy with flat RC). The basic problem is that there are very se- 
vere experimental bounds for the direct coupling of such a field with ordinary 
matter , and in the absence of such coupling the field will in general settle 
globally in the minimum of the potential leading to a homogeneous configu- 
ration that will not produce the desired effects. On the other hand, one could 
hope that, given the likelihood of existence of large black holes at the center 
of most galaxies, they would account for the nontrivial configuration of the 
scalar fields. Unfortunately these kind of situation is largely forbidden by the 



"Black Hole No hair theorems" for scalar fields [19-21|. These limitations 
severely reduce the types of models one can consider, in particular, there are, 
known to these authors, only three mechanisms that would account for the 
settlement of the scalar field in a non-trivial configuration in the absence of a 
direct coupling of the field with ordinary matter (or some other exotic matter 
which we will not consider because of the incremental number of hypothesis 
it involves): a) boson-star like clumps, b) spontaneous scalarization, and c) 
topological defects. Other models that lack these features have been consid- 
ered, for example in However, such models face two problems: first, they 
give rise to spherical configurations where the scalar field in consideration is 
singular "at the center" , and second, the resulting scalar field potential needed 
to account for the fiat RC depends explicitly on the value of the "tangential 
velocity" of stars at the flat region. That is to say, such a potential have to be 
adjusted differently for different galaxies. Needless to say that both problems 
clearly make those schemes unsuitable as models for the problem at hand. 



Concerning the case "a)" mentioned above, it has been analyzed in |23|. 
Their analysis focuses on cosmological and evolutionary considerations as well 
as the issues related to the conditions under which the assumption of long 
range coherency of the scalar field is justified, rather than the universal fea- 
tures of the galactic rotation curves. We will deal here with the other two 
cases "b)" and "c)". 
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The scenario "b)", namely the spontaneous scalarization (see Sec.V) 
is in some sense simpler because it involves a single scalar field in contrast of 
the various fields needed in the simplest versions of topological defects (e.g., 
global monopoles). Here the mechanism that allows for the nontrivial station- 
ary configuration of the scalar field is connected to a non-minimal coupling 
of the scalar field to the curvature. This results in the effective gravitational 
coupling becoming dependent on the scalar field. The point is that such a 
coupling allows for the reduction of the total energy of the configuration (in 
comparison with the corresponding configuration with the same baryon num- 
ber and no scalar field) for which the scalar field deviates from the trivial 
configuration by taking values that reduce the gravitational coupling in the 
regions of high matter density ||2^. Thus the model must incorporate from 
the onset the non-minimal coupling that seems to be needed to account at 
the same time for the correlations in the dark-luminous matter components 
(see |^4| and the discussion of the third scenario below). The disadvantage of 
this model, which is in fact shared by the first model (i.e., boson stars) p6| , 
is precisely the lack of resilience against black holes whose existence in most 
galaxies, if confirmed, would seem to preclude, through the no-hair theorems 



1 19 -21], the models based on this mechanism. 



The scenario "c)" is exemplified by the model of global monopoles |32| which 
have the notable feature of naturally leading to a 1/r^ energy density behav- 
ior which would naively account for the fiat rotation curves and which upon 
taking the symmetry braking scale to be the GUT scale would result in the 
correct order of magnitude for the galactic dark matter. Unfortunately, upon 
further examination of the simplest model severe problems arise, in particular 



the monopole configuration turns out to be repulsive |33|, and moreover the 
configuration would be too universal in the sense that it would be indepen- 
dent of the size of the galaxy thus defeating the hope for the correlation of 
dark to luminous matter over a range of galactic sizes. There is neverthe- 
less hope to overcome these problems by the consideration of slightly more 



complicated models |34|. In that work the simple monopole model was sup- 



plemented by the introduction of a non-minimal coupling between the scalar 



fields and curvature (see [34| and Sec. VI). This resulted in the restoration 
of gravitational attraction leading to regions of relatively fiat rotation curves 
and to the possibilities of the dark-luminous matter correlations arising from 
the fact that in these models the scalar potential V{^°'^a) (where stands 
for a triplet of scalar fields that characterize the global monopole) is replaced 
with the effective scalar potential V{^°'^a) + F{^"'^a, R) (here R stands for 
the scalar curvature of the space-time metric) whose minima depend on the 
amount of matter present trough the effect of the latter on R. The global 
monopole model has the additional advantage of resilience against the for- 
mation of black holes in the galactic centers, since their topological charge 
makes them immune to the devastating limitations imposed by the no hair 
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theorems. 



Despite the promising features of the model "c)", our intention in the present 
work is to take a step backwards and look at the problem from a more general 
point of view before embarking in the methodical study of a particular type 
of model. 

The article is organized as follows: in Section II, we analyze the generic 
form of the rotation curves of galaxies in a general relativistic context. In 
section III, we comment on the Newtonian approximation and on the embed- 
ding of the galaxy in the large scale space-time. In Section IV, we discuss 
the additional information that can be obtained about the metric from other 
considerations, specifically the deflection of light by the galaxy. Section V 
reviews the spontaneous scalarization scenarios. In Section VI, we review the 
non-minimally coupled global monopole model and discuss its shortcomings. 
Finally, in Section VII we offer a discussion and analyze the directions for 
further developments. 



II. ROTATION CURVES OF GALAXIES AND FREQUENCY SHIFTS 

The rotation curves (RC) provides the most direct method of analyzing 
the gravitational field inside a spiral galaxy. RC have been determined for a 
great amount of spiral galaxies | p^ , p^ . They are obtained by measuring the 
frequency shifts of light emitted from stars and from the 21 cm radiation from 
neutral gas clouds. 

In fact, since (apart from the central regions) the "tangential velocity" 
of rotation v remains approximately constant up to distances far beyond the 
luminous radius of these galaxies, a naive Newtonian analysis leads to the 
conclusion that the energy density decreases with the distance as and 
therefore that the mass of galaxies increases as m(r) « r. On the other hand, 
one could naturally question whether these large mass ought not to result into 
an important gravitational redshift. We will carry our analysis in a general 
relativistic setting and will see in the following sections that with standard 
assumptions about the matter content of the galaxy, the behavior of the RC 
indicates that the spacetime is not in general described by the standard form 
ds^ = -(1 + 2(^)dt^ + (1 + 2(^y^dr'^ + r'^dVt as can be initially thought (see 
for example [p^). 

In order to analyze the problem we will focus directly on what it is observed 
because only then will we be able to discuss models that do not let themselves 
to a Newtonian based inferences. This is an important point since the lack of 



understanding of it leads to erroneous conclusions |35| 



The observations of stars and gas in spiral galaxies show a shift Ztot in their 
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intrinsic spectra which includes the contributions of: 1) the cosmological ex- 
pansion (recession of galaxies), 2) the peculiar motion of the galaxy, 3) the 
thermal motion of atoms within the stars and gas, 4) the gravitational field 
within the galaxy and stars, and finally 5) the motion of the stars around the 
galactic center. 

When the "contaminating" effects from l)-3) are subtracted from the data, 
usually astronomers report the resulting z in terms of a velocity field v. Nev- 
ertheless it is instructive to make the analysis in terms of the quantities that 
are most directly observable: the z's. We perform this in order to keep track 
of the effect of the underlying assumptions, and to enable us to carry the 
analysis when these are not longer valid, as will be the case in some models 
we will consider. 

The starting point is to assume that stars behave like test particles which 
follow geodesies of a static and spherically symmetric space-time associated 
with sources that we do not specify for the moment. The most general line 
element of the space-time in these circumstances takes the form: 

ds^ = -N'^{r)df + A^{r)dr^ + r'^dO'^ + sin^ Odif^ . 

Next we consider two observers Oe and Od with four velocities 
spectively. Observer Oe corresponds to the light emitter (i.e., to the stars 
placed at a point Pe of space-time), and Od represents the detector at point 
Pd located far from the emitter and that can be idealized to correspond to 
"spatial infinity". 

Without loss of generality, we can assume that the stars move on the 
galactic plane 9 = 7r/2, so that n^; = {t,r,0,ip)E, where the dot stands for 
derivation with respect to the proper time of the particle. 

On the other hand, we suppose that the detector is static (i.e, Od's 4- velocity 
is tangent to the static Killing field ^), and so with respect to the above 
coordinates its 4- velocity is Uj^ = (f,0, 0,0)d- 

As usual, the consideration of the norm of the four velocity (u^u^ = — 1), 
gives, 

- 1 = -N\r)ii^) + A\r^) + r\^^) , 

The energy and the angular momentum per unit of mass at rest of the 
test particle are conserved quantities and can be written as 

E = -g^uE^'u' = N\i) , L = gf.u^P^u" = r\ip) , 

where e^, iJj^ denote the timelike and rotational killing fields of the metric (|l]) 
respectively. Introducing these constants of motion in the line element (|2|), 
we obtain 
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+ 1 



This equation shows that the radial motion of a geodesic is the same as 
that of a particle with position dependent mass and with energy m 
ordinary non-relativistic mechanics moving in the effective potential 
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As we mentioned, the RC of spiral galaxies are inferred from the red and blue 
shifts of the emitted radiation by stars moving in "circular orbits" on both 
sides of the central region |]l^ . The light signal travels on null geodesies with 
tangent k^. We may restrict k^^ to lie also in the "equatorial plane" 6 = it/ 2, 
and evaluate the frequency shift for a light signal emitted from O^; in circular 
orbit and detected hy Od- The conditions for circular orbits drVefj = and 
r = lead to 

r^drN/N 
1 - rdr-N/N ' 

1-rdrN/N ' 




The frequency shift associated to the emission and detection is given by, 

z = l , 

where 

= -k,,u'^\pc , 

the index C refers to emission or detection at the corresponding space-time 
point. 



Two frequency shifts corresponding to a maximum and minimum values are 
associated with light propagation in the same and opposite direction of motion 
of the emitter respectively (i.e., k^ = k^ = 0). Such shifts are the frequency 
shifts of a receding or approaching star respectively. Using the constancy 
along the geodesic of the product of the Killing field ^ with a geodesic tan- 
gent together with (^) and (|8|), and expressions (0) and (^), we find the two 
shifts to be, 

Nd {lT[rdrN{r)/N{r)f^) 
~ N{r) [1 _ rdrN{r)/N{r)]^/^ ' 

where N{r) represents the value of the metric potential at the radius of emis- 
sion r, and Nj^ the corresponding value of N(r) at r — > oo where the detector 



is supposed to lie. For asymptotically flat space-times No 1. However, for 
space-times generated by global- monopoles Nd — > (1 — a)^/^ (see the section 
VI). 

It is worth noting that in terms of the tetrads e(o) = ^^^St^^i^) ~ 
A-'l,e^2) = r-'^,e(3) = {rsm9)-%, the frequency shifts take the form, 

^± = i-^(iT^)r , 

where v := (Z]i=i,2,3('"{i)/'"(o))^)^''^ ^^^^ ^(a*) stands for the components of the 
star's four velocity along the tetrad (i.e., the velocity measured by an Eulerian 
observer whose world line is tangent to the static Killing field) and T = 

(l_^2)-l/2 

the usual Lorentz factor. Clearly, in the present case of circular 
orbits on the plane 9 = 7r/2, it results that v = n(3)/?X(o) = rdrN{r)/N{r). 
It is convenient to define the quantities: zn = \{z+ — Z-) andz^ = ^(z^ + z.) 
which are easily connected to the observations. From the expression (|lO|) we 
obtain, 

, , Nd [rdrN{r)/N{r)f^ 
ZD{r) - 



N{r) [l-rdrN{r)/N[r)f^ 



N{r) [1 _ rdrN{r)/N{r)Y^^ ' 

We note, for example, that {za — 1)^ — z]^ = {N{r) /No)^"^, and thus we could 
in principle recover N{r) directly from the observations. Then we can use 
this N{r) to recalculate za and Z£) from the above expressions and compare 
them with the measured values. This would be a test of the assumption that 
the dynamics is determined by the geodesies of a stationary metric, quite 
independently of the assumption of the dynamics of the geometry itself or for 
of the nature of the dark matter. 



III. GRAVITATIONAL FIELD IN THE DARK MATTER ZONE 

In this section we will use the form of the RC to obtain the spacetime 
metric and information on the matter content. The energy- momentum tensor 
must be diagonal and spherically symmetric, as dictated by the symmetries 
of the space-time (||), so we define: 

and r/;* = for /u 7^ v. 

We will for convenience introduce the following alternative form of the 
metric variables: 
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1 



2m(r) 



iV2(r) = exp[2zy(r)] 



Einstein's equations then read 

m{r) f 2m{r) 



dr 



1 



1 + 



mir) 



The equation of hydrostatic equihbrium resulting from the conservation of 
the energy-momentum tensor V^T^"^ = becomes 

drPr = -{drV)p + —{Pr " Pe) • 

V p J rp 

We note that the observations in spiral galaxies [17|, yield zd = v ^ 
constant and zjj >> za- Prom these conditions and from ( p^ ) we obtain 

1 1 

1 . 



and 



N{r) [l-rdrN{r)/N{r)]^^^ 



the value of v roughly ranges from 10 to 10 ^ depending of a particular 
spiral galaxy. The integration of ( [l9| ) gives, 

/ \ ''^ 

iV(r) 



r 
rn 



where rg is constant. 



Note that using (|2l|) in Eqs. (|17D and (|18D , we obtain a system of three 
equations [i.e., Eqs.(|l6|) — (|^)] for four unknowns (i.e, m, p, Pr, Pe). In the 
case of a perfect fluid, however, the four unknowns are reduced to three since 
Pr = Pe = P- This therefore constrains the equation of state p = p{p). On the 
other hand, for the case where the matter content is associated for example, 
with a scalar field, then p, Pr and Pq are not independent but are given in 
terms of the gradients of the field. This constrains the form of the scalar 
potential. 

We will look for a solution that satisfies the Newtonian conditions, 



Pr,e « P , 
r 

rn{r) « - , 

Airr^Pr « ni(r) . 



Under these conditions the Eqs. (p^) — (|T^) reduce to 
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dr'm{r) = Anr'^p 



m{r) 
-{drv)p 



using the expression (p|) for in (|2 
system, we obtain 



(p7|) with V = ln[N], and solving the 



m{r) 



p{r) 



Pr 



The solution corresponds to the relation ~ v'^p/2 which looks somewhat 
peculiar. If we view this as the equation of state of a perfect fluid, in the 
case of an ideal gas we would conclude that its temperature T is constant and 
proportional to v"^. The interpretation is that the dark matter represented 
by a perfect fluid is made of particles (ideal gas) that interact among them- 
selves strongly enough to maintain thermal equilibrium but do not interact 
in the same way with ordinary matter or with photons. The idea is then that 
the dark matter temperature determines its density profile and the spacetime 
metric, and through this, the rotation curves of the stars in the galaxy. One of 
the problems of this type of models is the need to explain why the isothermal 
configuration of the fluid does not extend to the center of the galaxy. 
Using Eqs. (^Tj) and ( [2^ ) in (|l|) we obtain the final result for the metric in 
this limit: 



2v^ 



dt^ + (1 



a 



^-1 



dr"^ + r'^{d9'^ + sin^ 



where a = 2v'^. We emphasize that this solution is only valid in the flat RC 
zone. We know that eventually this must be matched to a Robertson- Walker 
metric describing the universe or alternative we might use the asymptotically 
flat idealization for regions very far from the galaxy in question. To do so 
we must consider (p9|) as describing the spacetime geometry for r < Rq and 
the Schwarzschild metric for r > Rq, where Rq is the radius where the flat 
rotation curves end. The advantage of this approximation is that far from the 
galaxy the spacetime is Minkowskian a fact that facilitates, for example, the 
analysis of the propagation of light signals. 

Matching continuously the two metrics at Rq allows the determination of the 
integration constants: 



ds' 



ds' 



/ ^ \ 2»;2 

(l-a)(^— j dt^ + {l-ay^dr'^ + r'^{d9'^ + s{u^ed^'^) r < Rq 
1 _ ^\ dr^ + r\de^ + sin2 Odip^) r > Rq 



10 



where M = aRo/2. 

Here we are taking the view that the region of flat rotation curves, termi- 
nates in a narrow transition region where the behavior of the density changes 
from the to a constant that for simphcity we take to be zero. The point 
being that in the hmit in which the region is very narrow the metric will not 
change abruptly as we cross the region, but the derivative of the metric co- 
efficients will experience a jump. It would be interesting to consider various 
alternatives for the details of the interpolating regions. 

An alternative to the search of solutions satisfying the Newtonian conditions 
(PI) -(ID is to start from the Eqs.. (|ll)-(|l8|) without imposing the Newto- 



nian approximation but assuming that the dark matter is represented by a 
perfect fluid. Then again, the use of (|2TD with Pr = Pg allows the integration 



of (16)-dl|) as follows: 



-r , p{r) = — J— , Pr - 



2{l + 2v^-v^) 87r(l + 2z;2-w4)^2 ' Sir {I + 2v'^ - v^y"^ 

which results in an "equation of state" Pr = u^p/(4 — 2v'^). In practice 
« 1, so we recover the solution (p8[). 

At this point we can check whether the approximations we considered 
are self consistent. That is, we substitute ( pl| ) in the left hand side of (|19|) 
obtaining: 



[1 - rdrN{r)/N{r)]^'^ V " " ") 



The difference between this expression and 1 must be negligible in comparison 
to V which itself is of the order 10~^ — 10~^. This requires r/R^ to be neither 

2 

to large nor to small. To get an estimate we use the approximation X~'" ~ 
1 — v^\n{X) valid for |t;2ln(X)| << 1. These requirements are then: 



|ln(r/i?o)| « v-^ 

Thus the approximations are self consistent as long as —10^ << ln(r/i?o) << 
10^, which doesn't impose any practically relevant constraints for the case of 
the galaxies and the extent of the flat RC. 

It is worth noting that the form of the metric that we have obtained differs 
from what would be naively expected: 



ds"^ = - (1 + 2$) dt^ + (1 + 2$)-^ dr^ + r2((i02 ^ 



with $ representing the Newtonian potential. This form is often implicitly 
assumed (see and the fact that it is not appropriate for the region where 
the RC are flat lies at the core of the problems with the analysis of (see 



We also point out that if we assume that the flat RC extend indefinitely, 
the resulting spacetime is not asymptotically flat but rather is asymptotically 
flat but for a deficit angle (AFDA) |38|. In this context we stress that it is 



possible to consider such "idealized infinitely extended galaxies" as isolated 
objects in the framework of general relativity by replacing the asymptotically- 
flat framework by the framework AFDA |3^. One might want to embark in 
such considerations since in fact the RC remain flat to the farthest distances 
that can be observed. On the other hand there is a natural way to estimate 
an upper bound for the cutoff of such behavior. The idea is to consider the 
point at which the decaying density profile associated with the galaxy be- 
comes smaller than the average energy density of the universe. We call this 
bound R^ax • '^^^ value for -R^fax obtained by imposing the condition that 
the density at this point, provided by Eq. (^), is to coincide with the mean 
density of the universe: 



V 



2 

Max ) 



where pjj is the mean density of the universe. Then we have 



R^ 



Max 



Now we introduce the value pu in terms of the dimensionless Hubble param- 
eter h defined as Hq = 100/ikm/(secMpc) and Qu = pu/Pcvit 

PU = 1.8791 X lO'^^nuh^ kg/m^ = 2.78 x lO^^Quh'^MQ/Mpc^ . 

We obtain, 



<,, = 2.45xl06^O^^/'Kpc 



Taking h = 0.65 and v ~ (10 ^ — 10 ^), we have 

= 3.77 X (10^ - 102)O[;^/'Kpc 

Moreover, for 0,u = 1, 



^Max ~ 3.77 X (10^-10^) Kpc 
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On the other hand, the measured flat regions are about Rq ~ 2i?opt where 
i?opt is the radius encompassing 83 per cent of the total integrated light of 
the galaxy. We can take as a typical value Rq « 30Kpc, then Rq < R^^^- 



IV. BENDING OF LIGHT IN THE DARK MATTER ZONE 



One of the ways we could in principle explore the issue of whether the 
flat RC are the result of some form of unknown matter or the result of the 
change in the dynamical laws that governs the motion of particles, would be 
by studying the light deflection by galaxies. In particular, by studying the 
deflection of photons passing through the region where the RC are flat. Let us 
thus consider a photon approaching the spiral galaxy from far distances. We 
will compute the bending of light assuming the metric that has been matched 
with an asymptotically flat exterior, i.e., Eq. (l3C 

The bending of the light results 



A(f> = 2|(/)(ro) - (f)oc\ - vr 



(41) 



where (poo is the incident direction and vq is the coordinate radius of closest 
approach to the center of the galaxy. 



(ro)-<^oo= / A{r) 



-1/2 



dr 
r 



(42) 



The integral is split in two parts for the two domains of metric (E 



Ro 



A{r) 



-1/2 



dr . 
— + / Mr) 
r Jro 



ro 



-1/2 



The second integral is computed by expanding the integrand in powers of 
M/ro and M/r [| using (M) for r > i?o, 



Vo) N^r) 



1 + 2M ( - - — ) + 



^1 -1 



2Mr 



dr 
r 



(43) 



ro{r + ro] 



+ 



(44) 



and the results is 
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. fro\ M 

arcsm — — H 

\RoJ ro 



Ro 



1/2 



/ ^0 - ^0 
\Ro + ro 



1/2 



+ 



The first integral of (|4^) with (|36| ) for r < Rq gives 



A{r) 



r Y N^{ro) 



-1/2 



dr 



(1 - a)-V2 



(1-a)- 



-1/2 



(^;2 - 1) 



■ arctan 



2(l-u2) 



- 1 

2(1-1)2) 



-1/2 



dr 
r 

-1/2 



Finally, using (^) and ( 4f: ) in ( pD and then using (^) , the bending angle of 
light yields, 



, ro\ 2M 
2 arcsin ( — H 



2(l-a)-V2 
(?;2 - 1) 



arctan 



ro 



1/2 



-1/2 



/ -Rq - ro 
Vi?o + ro 



1/2 



where we took the limit 



lim^ 



arctan 



r\2{i-«2) 
ro/ 



-1/2 



vr 
2 



If we put ro = -Ro in (|47|), we obtain the standard result for the Schwarzschild 
metric with mass M = ai?o/2 and with Ac/) = (4M/ro) = 4 x 10^^. Fig.^ 
shows the bending angle of light A0 as a function of the parameter vq/Rq. 
If we take the impact parameter, ro, to be in the range of the measured flat 
regions by neutral Hydrogen measurements (HI): i?opt < ro < 2i?opt = Ro, 
then we have 1/2 < r^/ Rq < 1. In this case, the maximum value for Ac/) is ob- 
tained for the value tq/Rq = 1/2. Recently, the investigations for determining 
the radius of dark matter halos have gone beyond the HI measurements using 
satellite galaxies [27| or the weak lensing of background galaxies by foreground 



dark halos [28|. From these measurements, halo radii of more than 200 Kpc 
are inferred. For our galaxy 230 Kpc |29| and from satellite galaxies of a set 
of spiral galaxies show more than 400 Kpc [^0|. By taking Rq ~ 230 Kpc, we 
would have i?opt < ro < 15i?opt = Ro (where we have chosen i?opt = 15Kpc). 
In this case we would have 1/15 < ro/Ro < 1 and a value near to the max- 
imum in the fig.|l[ It would be interesting to explore the possibility to have 
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relevant observations in this context. 

We discuss next two of the three simplest scenarios in which the dark matter 
corresponds to coherent scalar fields. 



V. SPONTANEOUS SCALARIZATION 



As we mention in the introduction, the phenomenon of spontaneous scalar- 
ization in compact objects (notably in neutron stars) that arise in a class of 
scalar tensor theories of gravity is the other mechanism that allows the 
appearance of a non trivial scalar field in the absence of a direct coupling be- 
tween the scalar field and ordinary matter. The general feature of these kind 
of theories is a scalar field coupled non-minimally to gravity which leads to an 
effective gravitational coupling which depends explicitly on the scalar field. 
The non-trivial scalar field configuration appears when the object is compact 
enough so that the energy of the configuration for a fixed baryon number is 
minimized through a change in the value of the effective gravitational con- 
stant. That is, for a fixed baryon number, the energy of the configuration 
with a scalar field is lower than the corresponding configuration in absence of 
a scalar field |25|. An heuristic interpretation that is confirmed by the numer- 
ical results shows that, from a Newtonian point of view, the relevant quantity 
to be minimized is the combination GM instead of the total mass M. We 
observed that although the obvious additional contributions to GM are 
both positive and thus increasing the value of GM, their effect is more than 
compensated by the reduction of the value of the contribution GMbar ■, which 
is the leading term in GM. Thus there appears a nontrivial configuration 
of the scalar field which is associated with the minimization (at fixed total 
baryon number) of the the value GM p^]. 



Several problems arise if we want to use this mechanism to induce a non 
trivial configuration of a scalar field at the galactic scale. First, in the model 
studied so far we have seen that spontaneous scalarization occurs only if the 
object is compact enough, that is, if GM/R ~ 1/2 and needless to say that 
the galaxy as a whole does not satisfy this criteria (except perhaps at the 
center). If we assume that a large dense object lies at the center of the galaxy 
one would need some very unusual equation of state to overcome the stan- 
dard limits on the mass of these objects associated with the requirement of 
stability against collapse. But even if we were to assume such an object, the 
scalar field associated with the phenomena of spontaneous scalarization falls 
as 1/r (at least in the models considered so far) so it would not be relevant at 
the distances associated with the flat RC that lie at a distance of the order of 
kiloparsecs from the galactic center. Finally, the energy of the configurations 
with nontrivial configuration according to this phenomenon is smaller than 
that of the corresponding configuration in which the scalar field vanishes, thus 
the phenomenon seems to take us in the opposite direction as compared to 
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what seems to be required to explain the additional attractive effect on the 
test stars in the galaxy. If we wanted to consider extended objects other that 
neutron stars, it is not even clear how to build a sufficiently dense object. 
The only possibility would seem to be boson stars which now would act 
only as triggers of the spontaneous scalarization. These models would require 
to hypothesize two scalar fields, one providing the oscillating boson field of 
the boson star, and a second one providing the mechanism for spontaneous 
scalarization. On the other hand we must point out that although boson 
star masses are usually very small, when one introduces self-interactions their 
mass can be as large as 1{P''^ X^^'^Mq (for a scalar field mass m ~ 10~^ eV and 
a sufficiently large self-interaction constant A) Still we face the problem 
associated with the rapid fall off of the energy density associated with the 
scalar field, which would go like 



The hope here would be to consider alternative forms of the nonminimal cou- 
pling, with the possible introduction of various forms of self interaction terms 
for the scalar field, that would not only lead to spontaneous scalarization but 
to a rather different fall off behavior of the scalar field. Nevertheless, as was 
already mentioned, there is one very serious problem remaining with this type 
of scenario, and it is the issue of black holes. There is at the present time 
mounting evidence that there is at the center of most galaxies a very massive 
black hole, and in view of the no hair theorems for scalar fields [|l9|-21|, it 
seems clear that the phenomena of spontaneous scalarization does not have 
an analogy when the compact objects are replaced by black holes. Thus in 
those galaxies the scalar field would relax to the trivial configuration and thus 
any explanation of the RC based on that phenomena would cease to be oper- 
ative. There are, however, some small loopholes remaining in the black hole 
uniqueness theorems for the case of nonminimally coupled scalar fields which 
leave a ray of hope in this general direction, and which are currently under 
investigation E^. 



For the case of massless bosons (massless complex scalar fields), a New- 
tonian analysis leads to flat RC Unfortunately, in that work, the author 
neglect to note that the RC are not directly observable but only inferred from 
the corresponding light shifts. As it turns out, in that model the "gravita- 
tional" red-shift would be very large to the point that by ignoring it, the 
author is ignoring effects of the same order of magnitude as the ones that are 
being considered. Moreover, the law of composition of velocities used there 
to reproduce the RC is not valid. 

In the following section we analyze the case for the matter represented by 
global monopoles non-minimally coupled to gravity. 
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VI. GLOBAL MONOPOLES 



We will now consider one example of what we feel is at this time the most 
promising class of models: nonminimally coupled global monopoles. The main 
results of this section have been reported in 



Particle physics models predict the formation of topological defects during 
phase transitions in the early universe. The mechanism argued for the for- 
mation of these is the spontaneous breaking of symmetry of the model un- 
der consideration leading to a manifold of degenerate vacua with nontrivial 
topology. Topological defects can be classified according to the topology of 
the vacuum manifold. If the manifold of equivalent vacua, A4, contains un- 
shrinkable surfaces, TT2{Ai) 7^ /, then monopoles are formed. These can be 
classified into local and global monopoles depending on whether the symme- 
try broken is local or global. In the first case (gauge monopoles) the monopole 
configuration has finite energy concentrated in a small core and produces an 
asymptotically flat spacetime, while in the second case, the global-monopole 
configuration has a linearly divergent energy due to the long range Nambu- 
Goldstone field with energy density decreasing with the distance as r~^. As 
we have mentioned, this behavior is very appealing in view of the fact that 
this is precisely what seems to be required in a naive picture to provide a 
natural explanation for the flatness of the RC. 



It was shown by Barriola and Vilenkin that this linearly divergent "mass" 
has, at large distances, an effect analogous to that of a deficit solid angle a 
plus that of a tiny mass associated to the core of the monopole. Then, assum- 
ing the existence of a global monopole in a typical galaxy the total Newtonian 
mass contribution of the portion of the global monopole contained within rg^i 
(with Tgai ~ 15 Kpc) is found to be M ~ argai/2 « 10^^ GeV, where we took 
a typical grand unified value rj ~ 10^^ GeV, and where a = SttGt]^. This esti- 
mate turns out to be 10 times the total mass due to the contribution of 10^^ 
solar-mass in a typical galaxy (i.e., Mstars ~ 10^^ GeV). These numbers are 
again what is needed to account for the observations. Finally, if we assume 
that the field of the monopole extends on average a distance of ten galactic 
radii from the galaxy where the configuration presumably coincides with that 
of the monopole centered in the neighbouring galaxy, then M ~ lO'^*' GeV, 
which is 100 times that of the galaxy. This value leads to a contribution of 
the monopole to the total average density in the universe, which is of the or- 
der of magnitude predicted by the standard inflationary scenarios. Actually, 
it is the reversed argument that helps to place upper bounds on the density 
number of global monopoles present in the universe |3^. On the other hand, 
Harari and Lousto |33|, showed that the small effective mass ?n.core ~ 0.8a is 



in fact negative and produces a repulsive potential. They studied the motion 
of test particles in the spacetime of a global monopole concluding that there 
are no bound orbits. This result led thus to the unavoidable conclusion that 
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minimally coupled global monopoles are not good candidates to explain the 
RC despite the suggestive numbers and features considered above. 



Another problem is the fact that the monopole configuration is rather unique, 
in the sense that it is basically independent of the ordinary matter content in 
the corresponding galaxy, which conflicts with the fact that there is a rather 
large range of galactic masses for which the dark matter component is about 
ten times more massive than the ordinary matter component |39|. 

Recently, we have shown |34|, that coupling global monopoles non- minimally 
to gravity is possible to avoid the most unwanted features of the minimal 
case, notably, the lack of bound orbits, and the universality of the monopole 
configuration. 



Specifically we considered a theory of a triplet of scalar fields (j)^, a = 1,2,3, 
non- minimally coupled (NMC) to gravity with global 0(3) symmetry which 
is broken spontaneously to U(l). The simplest model of this kind is described 
by the Lagrangian 

where V{4>°'(j)a) is the usual Mexican hat potential V{4>°-4>a) = ■|((/>'^(/>a — v'^)'^- 
Equation (p9|) shows that the introduction of the coupling term is equiv- 
alent to consider an effective potential 

which explicitly depends on the matter content through R, and thus the con- 
tent of ordinary matter of the galaxy affects the location of the minima. This 
feature can thus help to avoid the scenario where the monopole configuration 
is universal, and opens the possibility to recover the correlation between the 
masses in the dark and ordinary matter components of the galaxy. 



167r 



R + F{R, 



In the following we show in detail how the non-minimal coupling leads to 
the existence of bound orbits. We will focus on the case where F{R, 4>'^4>a) = 
{^(p°'(l)a)R, where ^ is the NMC constant. The gravitational field equations 
following from the Lagrangian ( ^9|) can be written as 



where 



2^ 



rpflU 



Go 
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The equation of motion for the scalar fields is, 

We will only consider a metric describing spherical and static space-times (|l]) 
and study solutions of the gravitational and scalar fields equations describing 
global monopole configurations and the resulting space-time. Owing to the 
complexity of the resulting equations, we will perform a numerical analysis in 
terms of the following variables: 

u{r) = ln[iV(r)] , 
D{r) = v{r) — z^(0) , 



A{r)=[l-a- 



where 



The relevant Einstein equations take then the following form 
dm . 9 „ a 



dr 2Go 
du 2 f Gom a 



where 



E = N^Tl^ , 
is the effective total energy density. 

On the other hand, the Klein-Gordon equation can be written directly in 
terms of the energy momentum of the scalar fields: 



where 

S = ^cff i 1 

is the trace of the "spatial part" of T^^, which plays the role of an effective 
pressure. 



In the coordinates (|i|) this equation reads 
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r dr 
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+ 1^1 -a 
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2Gom 
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dr 



1 — a 



2Gom 



-1 



cos 9 8(1)'' 1 



The ansatz for a monopole configuration is 



r = vf{r)-, 

r 



with x"'x°' = r^, so that a monopole solution is found if / — > 1 at spatial 
infinity (i.e., r])- 



It is clear from the Eqs. (^^54), that the intermediary variables E and S [see 
Eqs.(|62[), (p^)] involve second order derivatives of the scalar field. However, 
we can eliminate such a terms from the gravitational field equations with the 
help of the ansatz for the monopole field and of the Klein-Gordon equation, 
and obtain "sources" containing at most first order derivatives of the scalar 
field. We also introduce the following dimensionless quantities 



f := r ■ T]}^!'^ , 
rh := m ■ GorjX ' , 



P-2- 



V 

A 

a 

Gcs 



V 



r 
A 

1 + 2^ A ' 
1 

1 + 2^Ap 



Go_ 



then the final form of the equations to be analyzed numerically is 



Airf^E , 
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+ ( l-a- — 
r 



(79) 



where 



E-S 



E 



AGeff 



87r(l + 24Ae2/2Geff) 
87r(l + 2AAe pG,s) 



2f 



-A2(9rfr + (1 + 120 + ir - ly + mrir - 1) 

'Am^ (l + 24AeV^Geff) + AU\f - 1) 

+ (^(^^/)' + (i + 8e + sAeV'Geff 



1 - s/^efG,^ 



Here E and 5" are dimensionless as in Eq. 

We note now that the sources of the differential equations contain only 
first order derivatives of the field variables and are thus suitable for numerical 
integration with a Runge-Kutta algorithm. 



(80) 



(81) 



A. Asymptotic expansions and boundary conditions 



Let us discuss the asymptotic behavior of global monopoles at the origin 
and spatial infinity in order to find the boundary conditions for the numerical 
integration. The regularity condition at r = on the metric requires 

?fi(0) = , dfm{{)) = -a/2 . (82) 

The boundary condition on P(r) is by definition 

f>(0) = . (83) 
The boundary condition on the scalar field at r = is (false vacuum), 

/(O) = . (84) 
Then one finds the expansions of the functions m(f), /(f) and P(f), at r = 0, 



m(f) 



a A 

—f H 

2 12 



1 



+ 3(8^ + l)fl 



+ 0{r 



f{r) = f,f + 0(f 



A 



(85) 
(86) 
(87) 
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where fc is determinated by the boundary conditions at spatial infinity. 
Let define M = m(oo). As we mentioned, we consider monopole configu- 
rations (i.e., / ^ 1 at spatial infinity: true vacuum). Then we have the 
asymptotic expansions at spatial infinity. 



m(r) = M + 



2fA 



o? ^ 8^(1 - a) 



A a 



+ 0(f 



v{f) = z>(oo) - ^ + 0(f-2) 
(1 — a)r 



We will impose for the asymptotic behavior of the metric, the standard 
asymptotically-flat-but-for-a-deficit-angle a spacetime (S.A.F.D.A a) (see 
1 38 1) and therefore the boundary condition on the function z^(f) at spatial 



infinity is 



i^(oo) = ^ln(l - q) 



The integration of the equations is performed by specifying the regularity 
condition on the scalar field at r = 0, 

af/(0) = /e . 

The value /c cannot be arbitrary, but must be so that / satisfy the appropriate 
boundary conditions at spatial infinity. This is enforced by the use of a 
standard shooting method |41]. 



We can compute the solution by integrating the equations in one step, 
from r = to a radius which is chosen conveniently so that / = 1 with a 
certain degree of approximation. The physical lapse N{r^ = e'^ at r = is 
calculated at the end of the numerical integration from ( |57| ) and ( |9l| ) 

j^(0) = ^ ln(l - a) - Poo , 

where the value Poo is obtained from the numerical integration. This ensures 
that at spatial infinity we recover the S.A.F.D.A a spacetime. 

We note that the ADM mass of the configurations of global monopoles, 
(see [^8| for a rigorous definition of the ADM mass for the case of spacetimes 
with a deficit angle) can be easily computed from the integral 

MADMa = M = lim^^oo rn(r) = (^ATTr^E{r) - dr . 
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B. Results. 



In order to obtain static configurations, we must impose the condition 
< a < 1 for the deficit angle. The behavior of the deficit angle a when 
A ~ is, 

a~ A-2eA2 + 0(A3). 

Case ^ > 1/2. The range allowed for the breaking scale is A G [0, cc) and 
for the deficit angle a € [0, ^). This is consistent with the condition a < 1. 
Figures shows numerical solutions for the case ^ = 2 satisfying the 

required boundary conditions. 

Case ^ < 1 /2. The range permitted for the breaking scale is A G [0, 1/(1— 
2^a)] and for the deficit angle a G [0, 1] . In this case the static configurations 
cease to exist when A > 1/(1 — 2^a). 

Figures |2|— shows numerical solutions for the cases = 0.3 and C = 
(minimal coupling case) satisfying the required boundary conditions. 

In the minimal coupling case Harari & Lousto estimated analytically 
the mass of the monopole as M ~ — 2a/3, and the size of the core 5 ~ 2. 
Then numerically showed that M ~ —0.75a. That is, they showed that the 
ratio Mja is practically insensitive to a (see figs.^,^. 

In the nonminimal coupling case, for a given ^ the situation is similar in 
that the ratio Mja is practically insensitive to a for a < 0.01. Actually, in 
this range of small a, that ratio depends weakly on the value of S,. However, 
for a > 0.01 the ratio depends strongly on and a [cf. figs.^j6|. This can 
be seen by performing an analysis similar to the one of but assuming the 
following approximation 



/ = 

and for the function rfi(f) 



(l+2CA)f 



if f < (5 
^ if f > 6 



m 



M if f > 5 



matching continuously the function rh and its derivative at f = 5, we ob- 
tain. 



2 A 2 1 /n 

M = — -— = — a(l + 2^A)-i/2 

3(l + 2^A)3/2 3 ^ ^ ^ 



5 = — = 2(1 - 2ea)i/2 

(1 + 2^A)V2 ^ ^ 
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where we have used the relation A = jji^^^- For the case ^ = 0, we repro- 
duce for M and 6 the values estimated analytically by Harari & Lousto ||33|| . 
We can see that in the case ^ > 1/2, we have A G [0,oo) and a € [0, ^) and 
therefore M — > when a — > these is confirmed numerically for the value 
^ = 2[cf. fig.^. Now we match continuously the function f at r = 6 and we 



obtain 



^ :(1 + 2CA)^/2 ^ 2 _ 2^^^-i/2 ^ (^Qo^ 



3^/3 3^/3 

then when a — >■ ^ we have that the value of fc diverges, this is checked out 
numerically for the value ^ = 2 [cf. fig-0- 

C. Geodesic motion in the spacetime of a global monopole. 



(101) 



In order to analyze the geodesic motion of test particles in the space- 
time generated by a global monopole we consider Eq.(^). In this case 
Nd = (1 - a)^/2 and the remaining gravitational potentials are given nu- 
merically. Figure |^ shows the effective potential for ^ = —2 and for different 
values of the other parameters; here we note the existence of a potential well 
and a non trivial minimum and thus the existence of stable circular orbits. 

For ^ > the effective potential (||) does not exhibit maxima or minima. 
An heuristic analysis that helps to understand the appearance of extrema 
of Veff when ^ < and their absence when ,^ > 0, is the following. We 
have that, drVefj = 2N^{r){-L'^ /r^ + L'^drv/r'^ + drv). So for drVeff = 
and then to find extrema, notably minima, it is necessary that dri^ > (or 
equivalently drN > 0). This means that N{r) should have an increasing be- 
havior as a function of r at least in some region away from the origin. Eq. 
d?^ ) provides the sign for the slope of A^. Near the origin / ~ fcr and then 
drV ~ — (A/12)(l -I- 24^)r/^ + 0{r'^). For dj-v > 0, ^ must be negative enough 
so that the coefficient (1 -|- 24^) is negative [cf. FigQ. On the other hand, 
for > 0, dri' < and therefore Veff has no minima (no bound orbits) [cf. 
Figs.^, 10 1 . This heuristic argument is confirmed by the rigorous numerical 



analysis from which the critical value ^crit that allows the existence of bound 
orbits is found to be ^crit ~ —0.15. 

Moreover, for L = which corresponds to 14// = -/V^, a peculiar situ- 
ation occurs in the cases where Veff has extrema, notably a maxima (e.g., 
for ^ = —2). The maxima and minima of T4// will correspond to the locus 
of unstable and stable stationary points where test particles are static (i.e., 
where particles do not feel any gravitational field). As seen from fig.|9|, we 
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appreciate that at the origin r = 0, test particles can be at rest in stable equi- 
librium, while at rmax where V(.ff is maximum test particles can be in unstable 
static equilibrium. In other words, r^ax separates two regions: one attractive 
and other repulsive [cf. the N{r) profile from Fig.^]. This strange behavior 
contrast dramatically from the cases of "conventional" gravitational sources 
like stars, planets, etc. where test particles are always attracted towards the 
source and where there are "no trivial" points at which they can remain static. 
In the case of minimally coupled global monopoles test particles are always 
repelled. 

Since bound orbit exists in the spacetime generated by nonminimal global 
monopole with suitable ^, we can compute the corresponding shift from 
Eq.(^) using the numerical solutions for and compare with the RC of 
spiral galaxies. Figure 11 depicts 2;+ (dashed line), Z- (solid line) and zo 
(dash-dotted line) as functions of f for the cases a = 0.43 (left panel) and 
a = 0.125 (right panel) respectively. We note that even for this very simple 
model the figures that would correspond to the rotation curves contain a rel- 
atively "flat region" within the values of r corresponding to stable orbits (i.e. 
the behavior of zd near its maximum). 



From these figures it is interesting to note that that in principle is as- 
sociated with a blue-shift, does not always correspond to a blue-shift, since 
there is a value ri, that separates z^ of being positive (blue-shift) or neg- 
ative (red-shift). This is easy to understand since the value z+ arise from 
a competition between the gravitational red-shift and the kinematical blue- 
shift. For slow particles (i.e, particles orbiting at "small" r) moving in the 
same direction as the emitted light, the gravitational barrier dominates over 
the positive contribution of the kinematical effects, and thus the frequency of 
emitted light has an overall attenuation. At orbits with radius r;,, particles 
are fast enough for the kinematical blue-shift to cancel the gravitational red- 
shift resulting in 2+ = 0. For example, from (^) and for v « 1, it turns 
that z+ < if V < 1 - iV(r)/(l - q)^/^. We note that in the non-minimal 
coupling case (1 — a)^/^ > N{r) in the regions where bound orbits exist. For 
larger values of r, z^ reaches a maximum and then starts decreasing since the 
gravitational barrier becomes larger (see Figj^) while the tangential velocity 
V [v = {rdrN/Ny^'^)] becomes smaller until reaching zero at the radius rmax 
where is maximum (see Fig.p^). 

Concerning Z-, this quantity is generically negative (i.e., it corresponds to a 
true red-shift). For instance, when v « 1, then 2_ ~ 1 — (1 — a)^/^(l + w)/A^ 
so z_ < if (1 — ay^^{l + v)/N > 1. This condition holds in most of the 
region of bound orbits since then iV < (1 - a)i/2 and f / (see Figs.!). 
However, moving away from the origin t> — > and grows to a maximum 
value where v = and 2;_ = 2;+ = 1 — (1 — a)^^"^ /N which can be positive. 
For the cases (values of ^) giving rise to bound orbits, it seems that always 
has a global maximum and then A'max > (1 — a)^/^. Therefore z± > at 
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Tmax- In fact the region of r where > is very narrow and corresponds 
to (1 - a)^/2(l + v)/7V < 1 [practically unseen at the scales of Fig.([TT|); this 
corresponds to orbits of small angular momentum]. 

To conclude this section, we mention that although the nonminimal global 
monopole model can repair the two main objections posed on the minimal 
model (namely, the bound orbit and the correlation between luminous and 
dark matter problems), there are still several improvements to perform in or- 
der that the quantitative predictions of this model fit reasonably well with 
the astrophysical data. Therefore, it is still very premature for any claim on 
this model as a realistic candidate for explaining the galactic dark matter and 
the corresponding rotation curves. The fact that the model is fully consistent 
in what regards the mathematical analysis (no singularities, no ad hoc pre- 
scriptions for the "tangential velocities" or for the metric), in addition to the 
numerical coincidences mentioned at the beginning of the section, provides 
some hope for pursuing a much more detailed study along this direction. 



VII. CONCLUSIONS 

The galactic rotation curves continue to pose a challenge to present day 
physics as one would want to understand not only the nature of the dark mat- 
ter that is associated with them but also the reason behind their universality 
(i.e., why is it distributed within a galaxy in a way that leads to almost flat 
rotation curves ?, and why is the amount of dark matter present in a galaxy 
so well correlated with the luminous matter ? [ p!^j50[ |). 

Models based in ordinary physical objects could already be facing prob- 
lems (depending on the exact value of the Hubble constant [^]) in view of 
the bounds that big bang nucleosynthesis impose on the baryon content of 
the universe. 

Models based on particle physics are the most commonly considered (usu- 
ally within a Newtonian scheme) but they need to address the nature and the 
distribution problems separately, leading to a larger number of hypothesis and 
surprinsing coincidences [pi]]. 

In view of the recent cosmological measurements and the theories that 
have been put forward to explain them [E^], one is naturally lead to consider 



alternative models based in the introduction of long range coherent fields [23[. 
In this work, we have given a review of various types of approaches to these 
questions indicating in each case the problems and advantages. 

We have argued that so far the most promising and simple approach would 
involve global monopoles with some sort of nonminimal coupling to gravity. 
This remains for the future to establish how far can this sort of ideas be 
pushed towards the goal of making a realistic and compelling model for the 
dynamics and evolution of galaxies. In particular any such model must also 
be studied in the context of cosmological perturbations, large scale structure 
and the CMB. In this regard we should point out that the simplest models of 
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topological defects as seeds for structure formation seem to be incompatible 
with the acoustic peak in the CMB anisotropics detected by Boomerang and 
Maxima |51]. However, all these studies have considered the simplest min- 
imally coupled models and it is unclear how would the models of the type 
being analyzed here behave in this respect. Finally, we should mention that 
the currently favored cosmological scenarios require at least two hypothetical 
components: the Cold Dark Matter (usually in the form of WIMPS) necessary 
for the structure growth and the dark matter in galaxies and clusters, and 
the cosmological constant A which provides the closure density (as required 
by inflation) as well as the repulsive component that seems to be required 
in order to account for the observations of the luminosity-distance of high 
red shift (type la) supernovae |52|. The fact that the non-minimally coupled 
monopoles exhibit both an attractive regime at short distances and a repulsive 
regime at large distances leads us to speculate whether these type of models 
can be used to explain the two aspects of the unobserved energy content of 
the universe in terms of a single hypothetical component. Needless is to say 
that all these aspects will require intense further exploration, which we hope 
to undertake in the near future. 



Acknowledgments 

U.N. is supported by a CONACyT postdoctoral fellowship grant 990490; 
M.S. and D.S. acknowledge partial support from DGAPA-UNAM Project No. 
IN121298 and from CONACyT Projects 32551-E and 32272-E. Authors thank 
the supercomputing department of DGSCA-UNAM. 



27 



REFERENCES 



V.C. Rubin, Science 220, 1339 (1983) 
M. S. Turner, Phys. Scr. T36, 167 (1991). 

J. R. Primack, Proceedings of the International School of Physics "Enrico Fermi", vol. 92, ed. 
N. Cabibbo. Bologna: Italian Physical Society, (1987). 
J. Wess, and B. Zumino, Nucl. Phys. B70, 39 (1974). 

E. W. Kolb, and M. S. Turner, The Early Universe, Addison Wesley, Redwood City, California, 
1990 

B. J. Carr, J. R. Bond, W. D. Arnett, Astrophys. J. 277, 445 (1984). 

C. Alcock, et al., Astrophys. J. 486, 697 (1997). 

P. J. E. Peebles, Principles of Physical Cosmology, Princeton University Press, Princeton, 
(1993). 

Steven Weinberg, Gravitation and Cosmology: Priciples and Applications of the Ceneral The- 
ory of Relativity (Jonh Wiley and Sons, Inc, 1972). 

C. J. Copi, D. N. Schramm and B. M. Tinsley, Science 267, 192 (1995). 
J. A. Sellwood,and A. Kosowsky, |astro-ph/0009074 
A. Finzi, Mon. Not. Roy. Astron. Soc. 127, 21 (1963). 
R. H. Sanders, Astron. Astrophys. 136, L21 (1984). 
M. Milgrom, Astrophys J. 270, 365 (1983). 

n 

p. D. Mannheim, Astrophys J. 479, 659 (1997); ibid, ApJ, 419, 150 (1993). 
M. Persic, P. Salucci, F. Stel, M.N.R.A.S. 281, 27 1996. 

V. C. Rubin, W. K. Ford, and N. Thonnard, Astrophys. J. 238, 471 (1980); V. C. Rubin, W. 
K. Ford, N. Thonnard, and D. Burstein, ibid. 261, 439 (1982); 289, 81 (1985). 
V. B. Braginskii, V. I. Panov, Zh. Eksp. Teor. Fiz. 61, 873 (1971). 
M. Heusler, J. Math. Phys. 33, 3497 (1992); ibid Class. Quant. Grav. 12, 779 (1995). 

D. Sudarsky, Class. Quant. Grav. 12, 579 (1995). 
J. D. Bekenstein, Phys. Rev. D 51, R6608 (1995). 

T. Matos, F. S. Guzman, and Dan'o Nuiiez, Phys. Rev. D 62, 061301 (2000) 
P. J. E. Peebles, |astro-ph/0002495| ; A. Riotto and I. Tkachev, |astro-ph/0003388| ; J. Goodman 
|astro-ph/0003018| 

T. Damour and G. Esposito-Farese, Phys. Rev. Lett. 70, 2220 (1993); ibidFhys. Rev. D 54, 
1474 (1996); ibidFhys. Rev. D 58, 042001 (1998). 

M. Salgado, D. Sudarky and U. Nucamendi, Phys. Rev. D 58, 124003 (1998). 
I. Peha and D. Sudarsky, Class. Quant. Grav. 14, 3131 (1997). 

D. Zaritzky, R. Smith, C. S. Frenk and S. D. M. White, Astrophys. J. 405, 464 (1993); D. 
Zaritzky and S. D. M. White, ibid. 435, 599 (1994). 
T. Brainerd, R. Blanford and I. Small, Astrophys. J. 466, 623 (1996). 

A. S. Kulessa and D. Lynden-Bell, Mon. Not. R. Astr. Soc. 255, 105 (1992); C. S. Kochanek, 
Astrophys. J. 457, 228 (1996). 

D. Zaritzky, R. Smith, C. S. Frenk and S. D. M. White, preprint |astro-ph /96 1 1 199| . 
Those contributions are: that which is associated to the energy density of the scalar field 
and the other which is associated to the reduction of the negative binding energy due to the 
reduction of G due to a nonzero value of the scalar field: G ~ 1/(1 + 167r^(/)2) with ^ > 0. 
[32] M. Barriola and A. Vilenkin, Phys. Rev. Lett. 63, 341 (1989). 



28 



D. Harari and C. Lousto, Phys. Rev. D 42, 2626 (1990). 

U. Nucamendi, M. Salgado, and D. Sudarsky, Phys. Rev. Lett. 84, 3037 (2000). 
A. Edery, Phys. Rev. Lett. 83, 3990 (1999). 

J. D. Bekenstein, M. Milgrom, and R. H. Sanders, Phys. Rev. Lett. 85, 1346 (2000). 
W. A. Hiscock, Phys. Rev. Lett. 64, 344 (1990). 

U. Nucamendi and D. Sudarsky, Class. Quant. Grav. 14, 1309 (1997). 
R. M. Wald (private communication). 

D. P. Bennett and S. Hong Rhie, Phys. Rev. Lett. 65, 1709 (1990). 

W. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterling, Numerical Recipes, The art of 

scientific computing, Cambridge University Press, 1986, p. 578. 

M. Colpi, S. L. Shapiro, and L Wasserman, Phys. Rev. Lett. 57, 2485 (1986) 

W. L. Freedman et al.. Nature 371 757 (1994); M. J. Pierce et al.. Nature 371 385 (1994) 

J. Binney and S. Tremaine, Galactic Dynamics, Princeton University Press, Princeton, (1987). 

E. Seidel and W.-M. Suen, Phys. Rev. Lett. 72, 2516 (1994) 

F. E. Schunk, and E. Mielke, Proc. 8th M. Grossmann Meeting, (ed. T. Piran), World Scientific, 
p. 1441, Singapore, 1999. 

I. Peha and D. Sudarsky, In preparation. 

F. E. Schunk, in Proceeding of the 8th Marcel Grossman Meeting (ed. T. Piran), World Scien- 
tific, p. 1441, Singapore, 1999; F. E. Schunk, |astro-ph/9802258 



T. Barreiro, E. J. Copeland, and N. J. Nunes, Phys. Rev. D 61, 127301 (2000) 
R. B. Tully and J. R. Fisher, A&A, 54, 661 (1977). 

P. de Berdarnis et al. Nature 404, 955 (2000); S. Hanany et al, Atrophys. J 545, L5 (2000); 
U.L. Pen, U. Seljak, and N. Turok, Phys. Rev. Lett. 79, 1611 (1997); R. Durrer, M. Kunz, 
and A. Melchiorri, Phys. Rev. D 59, 123005 (1999). 
[52] A. G. Riess et al.. The Astronomical Journal 116, 1009 (1998); S. Perlmutter et al, Astrophys. 
J. 517, 565 (1999) 



29 



FIGURES 





logio[l + ^AJ 

FIG. 2. The figure depicts the global monopole field /(r) for ^ = 0, a = 0.795 (solid line), 
^ = 2, a = 0.1 (dashed line), and ^ = 0.3, a = 0.79 (dash-dotted line). Here Tc = (r/A^/^)-^ 
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FIG. 3. Mass profile m(r) for ^ = 0, a = 0.795 (solid line), ^ = 2, a = 0.1 (dashed line), and 
^ = 0.3, a = 0.79 (dash-dotted line). Asymptotically this quantity provides the ADM mass of the 
configuration. Here rric = (Go'?A^/^)~^. 
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FIG. 4. Metric potentials N{r)/{1 - a)^/^ (solid lines), A{r){l - a)^/^ (dashed lines), and the 
product AN (dash-dotted lines) for three different configurations. The left panel corresponds to 
^ = 0, a = 0.795, the middle panel to ^ = 2, a = 0.1, and the right panel to ^ = 0.3, a = 0.79. 
Note that ^A'^ aproaches to unit "outside" the monopole core, while the metric potentials N and 
A tend to the asimptotically-flat-but-for-a-deficit-angle values far from the origin. 
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FIG. 6. ADM mass-deficit angle rate of configurations for different values of the deficit angle 
for ^ = (solid line), ^ = — 2 (dashed line), ^ = 0.3 (dash-dotted line), and ^ = 2 (dotted line). 
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FIG. 7. Derivative of the global monopole at r = for different values of the deficit angle for 
^ = (solid line), ^ = —2 (dashed line), ^ = 0.3 (dash-dotted line), and ^ = 2 (dotted line). 




FIG. 8. Functional dependence of the effective potential V^ff vs r for the case ^ = —2, 
a = 0.43 (left panel) and a = 0.125 (right panel). For each configuration we show three values of 
the angular momentum: L = 4 (solid line), L = 1 (dashed line) and L = 0.3 (dash-dotted line). 
Here rc = (r/A^/^)"^. 
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FIG. 9. Same as figj| for ^ = -2, a = 0.43. 
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FIG. 10. Same as fig.| for ^ = 0, a = 0.795 (left panel) and for ^ = 2, a = 0.1 (right panel). 



36 





FIG. 11. Functional dependence of Zj^ (dashed line), z_ (solid line) and (dash-dotted line) 
vs f for circular orbits in the case ^ = —2, a = 0.43 (left panel) and a = 0.125 (right panel). The 
asterisk depicts the location of the radius beyond which the stable circular orbits cease to exist 
(r ~ 9rc). 
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FIG. 12. Tangential velocity v = {rdrN/NYf^ in units of c for ^ = -2, a = 0.125 (solid line) 
and a = 0.43 (dashed line) . The asterisk depicts the location of the radius beyond which the stable 
circular orbits cease to exist. 
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